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Abstract 

We use phase space method to study possible consequences of fixed points in flat FLRW 
models. One of these consequences is that a fluid with a finite sound speed, or a differentiable 
pressure, reaches a fixed point in an infinite time and has no finite-time singularities of types I, 
II and III described in hep-th/0501025. It is impossible for such a fluid to cross the phantom 
divide in a finite time. We show that a divergent dp/dH, or a speed of sound is necessary but 
not sufficient condition for phantom crossing. We use pressure properties, such as asymptotic 
behavior and fixed points, to qualitatively describe the entire behavior of a solution in flat 
FLRW models. We discuss FLRW models with bulk viscosity i] ~ p r , in particular, solutions 
for r = 1 and r = 1/4 cases, which can be expressed in terms of Lambert-W function. The last 
solution behaves either as a nonsingular phantom fluid or a unified dark fluid. Using causality 
and stability constraints, we show that the universe must end as a de Sitter space. Relaxing the 
stability constraint leads to a de Sitter universe, an empty universe, or a turnaround solution 
that reaches a maximum size, then recollapses. 
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1 Introduction 



Various cosmological observations [1, 3, 4, 5] have provided us with a strong evidence for 
accelerating expansion of the universe. Although, the component that causes this acceleration 
is not known yet, the best model that fits dark energy and other cosmological data is the ACDM 
model, which is a Friedmann-Lematre- Robertson- Walker (FLRW) universe with a cosmological 
constant. Unfortunately, this model does not provide a physical picture of dark energy. In 
order to describe dark energy in the realm of general relativity, we need to consider some exotic 
fluid with an unusual equation of state that has a negative pressure and violates the strong- 
energy condition (see for example [9, 10]). The existence of this exotic fluid not only opens 
the door for reexamining the constituents of our universe but also evades some nonsingularity 
theorems through relaxing the strong-energy condition. This revives interest in nonsingular 
cosmologies, especially those describing the universe in early and late times. Having a more 
general class of equations of state, that violates the strong-energy condition, creates a wider 
class of singularities beyond that of Big Bang and Big Crunch [6, 7, 8]. In FLRW models these 
singularities are classified as follows [25]: 

• Type I ("Big Rip"): t — > t s , a — > oo, p — > oo, and \P\ — > oo 

• Type II ("Sudden"): t — > t s , a — > a s , p — > p s , and \P\ — > oo 

• Type III : t — > t s , a — > a s , p — > oo, and \P\ — > oo 

• Type IV : t — > t s , a — > a s , p — > 0, and \P\ — > 0, but higher derivative of H diverges. 

It is of interest to build cosmological models, which are free from the above singularities. 

In the last decade, there have been several proposals to describe dark energy/matter or dark 
energy alone as a single barotropic fluid in FLRW models with various equations of state (EoS) 
such as; Chaplygin gas, Van der Waal, linear, and quadratic EoS[13, 14, 15, 16, 17, 18, 19, 20, 
21, 22, 23, 24]. Notice that, most of these models have fixed points, or de Sitter solutions. 

In this article we use a phase space method to study general solutions of single fluid FLRW 
models with fixed points and a pressure p{H), where H is the hubble parameter. We discuss 
possible consequences of having these fixed points. Some of these consequences are; (i) if we 
describe our universe as a single component fluid and model the late time acceleration by a 
future fixed point, then the resulting cosmology does not have future-time singularities of types 
I, II and III described in [25], (ii) cosmologies with a future and a past fixed points are free 
of types I, II and III singularities, (in) one can use a simple argument to show the phantom 
divide [27, 28], or in a single fluid FLRW models it is impossible for a physical solution to cross 
the phantom divide line in a finite time, and (iv) in these models, the only way to get bounce 
solutions is to have a nonvanishing pressure as p — > 0. The phase space method can be used to 
construct nonsingular late time model, in particular, unified dark fluid (UDF) and dark energy 
models. We use this qualitative analysis to describe the entire behavior of a FLRW cosmology 
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with bulk viscousity rj ~ p r , presenting two exact solutions, with r = 1 and r = 1/4 expressed 
in terms of Lambert-W function. The last solution describes either a nonsingular phantom 
dark energy or a unified dark fluid model. At the end of section 4 we list possible scenarios for 
the future of our universe. 



2 Flowing to a Fixed Point 

In this section, we use a phase space method in flat FLRW models with fixed points to argue 
that a fluid with a continuous and differentiable pressure always reaches a fixed point in an 
infinite time and has no finite-time singularities of types I, II and III described in [25]. 
Let us start with a FLRW universe, with an equation of state 

p=p(H), (1) 

where the pressure p(H) is a continuous function of the hubble parameter H. Using the unit 
convention, 8nG — c — 1, Einstein field equations lead to Friedmann equation and Raychaudhri 
equation 

H = -H 2 -l(p + 3p)) H = a/a. (3) 
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For energy-momentum conservation, we obtain 

P = SH(p + p) (4) 

which is not independent of Eqn. (3). In this work we are interested in the flat FLRW case, 
the spatially curved case will be discussed elsewhere [30]. For a flat FLRW universe, Eqn. (3) 
becomes 

H = -\(p(H)+p) = f(H), (5) 

which can be expressed in terms of a dimensionless hubble parameter and time h = H / H* and 
r = c H* t, as follows 

- = F(h). (6) 

Where H* is a parameter that depends on the equation of state parameters and c is some 
number. A solution, h{r) of Eqn. (6) is subject to an initial condition h(0) = ho. 

Motivated by the fact that many of the proposed models for dark energy and unified dark 
matter/energy models do have fixed points, we assume the existence of fixed points for Eqn. (6), 
which are the zeros of the function F(h) = —1/2 (p + p). Let us call the zeros of F(h), h±, h 2 , 
where hi < h 2 < h 3 ... 
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Fixed points are classified according to their stability as follows; stable, unstable, or half- 
stable, depending on the sign of their tangents as shown in Figure (l) 2 . Unstable fixed points 
are represented by arrows emanating out of them and stable points are represented by arrows 
pointing toward them as shown in Figure (1). Half-stable points are stable from one side and 
unstable from the other side or vice versa. 



The arrows determine how the solution develops with time. Notice that, a fixed point 
satisfies equation (6), i.e., h(r) = hi, therefore, it is a solution. This constant solution is 
nothing but a de Sitter space. If the system started from an initial value, ho = hi, i.e., at a 
fixed point, then it will remain at this point forever. But when the value of ho is close to that 
of hi, the solution h(r) develops towards hi, if hi is a stable fixed point, or away from it, if the 
point is unstable. This technique has been used in literature to study particular equations of 
state and their solutions (e.g., [18]), but here we try to keep our discussion general assuming a 
pressure p(H). 

Using this tool, one can determine how a solution behaves upon knowing the nature of these 
fixed points. For example, if we start with a value for ho between the stable and the unstable 
fixed points in Figure (1), then the solution will develop to the left, i.e., will take values with 
smaller h till it reaches the stable point. Therefore, the solution in this case interpolates between 
two de Sitter spaces. If ho takes values larger than the unstable fixed point in Figure(l), then 
the solution develops to the right and might reach a singularity if the time to reach h — > oo is 
finite. 

Fixed points and the asymptotic behavior of F(h) enable us to predict the behavior of the 
system without knowing the form of the solution. But, in order to have a reliable qualitative 

2 Hcre I am following the terminology of Strogatz [26]. 



Fixed Points 




Figure 1: Types of fixed points 
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description of a solution we have to know how long it takes to reach a fixed point or a singu- 
larity, i.e., a point where h — > oo. Using the argument below we will see that the continuity 
and differentiability of the pressure, p{h) determine if the time to reach a fixed point is finite 
or infinite. Also, in the coming section we use the asymptotic behavior of F{h) to determine 
the time taken by a solution to reach a point where F{h) is diverging. We are going to use 
pressure properties to describe the behavior of a solution qualitatively without the need for 
exact or approximate solutions. 

Here we argue that for a flat FLRW fluid with a pressure p(h), which satisfies; 

i) F{h) is continoues and differentiable, and 

ii) F(h) has a future fixed point hi, i.e., h < hi if F(h ) > 0, (or h > hi if F(h ) < 0), 
there exists a unique solution h{r) which is defined for times r > 0, and has no future finite- 
time singularities of types I, II and III. 

Notice that, the first assumption is needed to ensure the existence of a unique local solution 
around some initial value h by the existence-uniqueness theorem. In addition, it leaves the 
solution free from type-II singularities, which might contradict casuality since the sound speed 
dp/ dp = c 2 s will diverge as well. The second assumption ensures the extendibility of this local 
solution to all values r > while keeping the hubble parameter bounded. 

Let us start our argument by integrating equation (6), assuming an initial value h(0) = ho, 
then one obtains 

r h dy 



GW= Lwr T (7) 

Since F(ho) ^ 0, its either F(ho) > or F(h ) < 0, let us choose F(h) > 0, and ho < hi 
from ii). As a result G{h) is a monotone near h - For any solution 0(r), we have 

G(0(r)) = r. (8) 

Since G{h) is a monotone near ho, the above relation can be inverted 

G(r)-^^r), (9) 

where G -1 is the inverse map of G. This is a local solution by construction (unique since F(h) 
is differentiable by the existence-uniqueness theorem) that can be extended by looking for a 
maximal interval in which G(h) is a monotonic function. Since F(h ) > and remains positive 
for the values h < h < hi, then G(h) is a monotone in this interval. This imply that the 
maximum interval in which we can extend the solution is [ho, hi]. The solution can be defined 
for all values r > if 

ft i dh 

* W) = °°- (10) 
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Now, let us show that r + = oo, when i) and ii) are given. Since F(h) is differentiable, the slope 
of the tangent at any point is finite. Let us choose a number M < F'(h),Vh G [h ,hi\. The 
existence of M enables us to define a linear function Y(h) such that 

Y(h) > F(h), (11) 



hl dh [ h i dh 



where — M (h — hi). This leads to 

ho ~~ A F'(h min ) (h - hi) 

oo. (13) 



therefore, 

rfti dh 



Iho F(h) 

We have shown that for the first-order system of Eqn.(6) there exists a unique solution, h(r) 
defined for times r > if the above two assumptions are satisfied. Also, it takes the solution 
an infinite time to reach the future fixed point hi. It is clear that the solution is bounded, i.e., 
h(r) G [ho, hi] for times r > 0. Therefore, the density, p(r) is bounded for times r > 0, as a 
result, there is no future singularities of type I and III. The pressure p = —2F(h) — 3 h 2 H* 2 
is bounded too in this interval, since F(h) and h are bounded 3 , which means, no future singu- 
larities of type II in these spacetimes. As a result this class of FLRW solutions are free from 
future finite-time singularities of type I, II and III . 



The above argument can be generalized by relaxing the differentiability condition, in this 
case we have one of the following; 

a) F'{h) has a jump discontinuity: F(h) is not continuous, but its left and right derivatives 4 
are finite for h G [ho,h]\. In this case, one can still construct the linear function Y{h), which 
can be used to show that the time to reach hi is always infinite and the solution has no future 
finite-time singularities of types, I, II, and III. 

b) F(h) has infinite discontinuity: F(h) is not continuous and at least the left or right derivative 
of F(h) is infinite for some h G [ho, hi]. This case is not physical, since the divergence of F'(h) 
leads to a divergent speed of sound, dp/ dp = c 2 , which has to be less than unity for the model 
to be causal. We are going to discuss this issue in more details in the coming section when 
we discuss phantom crossing, but through out this paper we are going to assume no infinite 
discontinuities in F'{h). 

Since any continuous function on a closed interval is bounded 
4 Left derivative is denned as cL j^ a ' ) = lim /l ^ a - F ^_^ ( - a - ) and the right derivative is defined in a similar way. 
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3 Consequences of Fixed Points 



Here we discuss consequences of fixed points and how can we use this dynamical method to 
describe the entire behavior of a single fluid in flat FLRW without knowing the form of the 
solution. 

3.1 Direct Consequences 

Consequences of using phase space method to study fixed points can be listed as follows; 

i) If the late-time behavior of our universe is described by a single fluid, as in unified dark 
fluid models, and the late-time acceleration is developing towards a de Sitter universe, then 
there is no future-time singularity of types I, II and III in this solution. Since the pressure 
p(h), is a differentiable function of h, or at most has finite discontinuities, it takes the universe 
infinite time to reach the de Sitter space. 

ii) If we combine a future fixed point with a past fixed point, by applying the argument in 
section 2 we get a nonsingular solution, which is free from types I, II and III singularities. The 
time taken by the solution to reach hi or to come from h 2 starting from an initial value h is 
infinite. In addition, the hubble parameter h{r) interpolates between h 2 and hi in a monotonic 
manner. 

in) As one might notice, the solutions we have so far, still, admit a weaker type of singularity, 
namely, type-IV in the classification given in [25]. But even these weaker singularities can be 
avoided by requiring p(h) to be a smooth function, i.e., a C°° function. One can show that as 
follows; The n-time derivative of h(r) can be written as; 



which means, unless F{h) or one of its derivatives (up to the (n— l)th derivative) is divergent, 
h^> is always bounded in \h 2 , hi]. This leads to the conclusion: In a flat FLRW universe, if a) 
p(h) is a smooth function, i.e., arbitrarily differentiable, and b) F(h) has a future and a past 
fixed point, then the spacetime is free from singularities of types, I, II, III and IV [25] when 
h e [h 2 , hi]. 

iv) A single fluid in flat FLRW with a pressure p(h) admits bouncing solutions only if there 
are no fixed points between h — and h = hi, i.e., F(0) ^ 0. These solutions either have a 
bounce or turnaround at h = depending on sign of .F(O). At h — 0, a/a — F(0), therefore, if 
F(0) > it will be a bounce and for F(0) < it will be a turnaround. 

v) Another consequence of this phase space method is a simple and transparent way to show the 
no-go theorem of phantom divide [27, 28] (see also, [39, 29]), which can be stated as follows; In 




(14) 
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a single fluid FLRW cosmology it is impossible for a causal solution to go from a region where 
ui(h) < —1 to another with u>(h) > —1. In other words, a solution in one of the mentioned 
regions has no access to the other region. Let us explain this in more details using the analysis 
we have in section 2. First, let us assume that the pressure p(h) is differentiate. In this case, if 
a solution approaches a fixed point, where w(h) = —1, starting from a region where w(h) < —1, 
then, it will spend an infinite time to reach it, as a result it will never cross it. If p(h) is not 
differentiable, then p'(h) has either a finite or an infinite discontinuity. If the discontinuity is 
finite, the time to reach the crossing point or the fixed point is infinite, as we showed in the 
previous section, therefore, the crossing will not occur. If p'(h) has an infinite discontinuity, 
a solution will reach w(h) — — 1 in a finite time, but in this case the solution is not causal. 
Although, the solution is not causal, it is not clear if it is going to cross the phantom divide or 
not. To further investigate this case, let us list the conditions on F(h): 

i) lim h _ hl F(h) = 0, 

ii) lim/j^/jj dF(h)/dh = ±oo and 

iii) t = JlS l/F(h) < oo. 

A class of functions which satisfy the above conditions is F(h) = F (hi — h) s , where < s < 1 
and F Q > 0. A solution is given by the following expression; 

h(r) =h 1 - [(h - h o y- s - F (1 - s) t] ^ , r < r* 

= h T>T* (15) 

where, h(0) = h and r* = (hi — ho) 1 ^ 3 / F (1 — s). It is clear from Eqn. (15) that the solution 
stays at the fixed point for t > t*. Also, notice that hi is a stable fixed point. One can see 
that by assuming a small perturbation away from hi, i.e., h(r) — hi + S(t), it leads to 

S(t) = - [(h - h ) l ~ s - F (1 - s) t] ^ . (16) 

As a result, if we extend the definition of F(h) for h > hi, e.g., F(h) = —F (h — hi) s as shown 
in Figure (2)-a), we will not have a phantom crossing. But if we define F(h) = ±F (hi — h) s , 
i.e., a double valued function as in Figure (2)-b), the solution can cross the phantom divide 
in a finite time. Therefore, in addition to its infinite discontinuity F'(h), need to be a double 
valued function to have a phantom crossing solution. Our conclusion is that a general causal 
solution of any flat FLRW model, with a continuous pressure p(H), can not cross the phantom 
divide line in a finite time. It is clear that infinite discontinuity of F(h) is necessary but not 
sufficient for a phantom crossing. 

3.2 Describing a Solution Qualitatively 

As we mentioned earlier, to have a complete qualitative description of a general solution in 
flat FLRW cosmology we need to know the fixed points as well as the asymptotic behavior of 



Figure 2: a) F(h) = (1 - h) 1 ' 2 , for h < 1, and = -(/* - l) 1 / 2 , for h > 1. b) F(h) = ± (1 - h) 1 ' 2 , for 

h<l. 



F(h). The later property enables us to determine the time to reach a point where F(h) — >■ ±oo, 
starting from some initial value /i(0) = fro- Let us first show the relation between the asymptotic 
behaviors of F(h) and finite-time singularities. 

Asymptotic Behavior of F(h) and Singularities: The dimensionless Hubble parameter h, 
in a flat FLRW cosmology is controlled by a one-dimensional phase space evolution function 
F(h). In this dynamical system a solution develops towards either a fixed point, where F(h) — > 
0, or a point where F(h) — > ±oo. A solution approaching a point where — > ±oo does not 
necessarily mean that the it has a finite-time singularity. It is crucial to know how fast F(h) 
reaches infinity. This enables us to determine if the singularity is reached in a finite time or 
not [32]. Considering the integral 



it is easy to see that if lim^ioo F(h) ~ h s , where < s < 1, the integral diverges. As a result, 
the solution takes an infinite time to reach the singular point, therefore, it has no finite-time 
singularities. One can observe that if F(h) grows as a linear function or slower 5 , as h — > ±oo, 
the solution will have no finite-time singularities. One can show this rigourously following the 
same argument in section 2. An example of a singular asymptotic behavior is a quadratic 
function F ~ h 2 , which leads to a Big Bang singularity. 

It is intriguing to notice that there is a class of F(h) that grows faster than a linear function, 
but still, does not have finite-time singularities. An example of this is F(h) ~ h ln/i, or 
h In h In In h, and so on [31]. These functions grow faster than a linear function, but lead 
to time r(h), that depends on h logarithmically, therefore, leads to nonsingular solutions. In 
general, if F(h) can be expressed as F(h) = g/g', where g(h) is any function such that g(h) — > oo 

5 This is a known fact in the literature see for example [32] 




(17) 
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as h — ?■ ±00, then r(h) will diverge logarithmically as h — > ±00, which leads to a nonsingular 
solution. To conclude, except for F(h) with a special form (as in the above mentioned cases), 
the solution reaches a point where F(h) — > ±00 in a finite time if F(h) grows faster than a 
linear function. 

Qualitative Description: In this section we take the pressure p(h) to be a continuous function 
of h, or at most has finite discontinuities. Therefore, we are not going to allow any infinite 
discontinuities for p'(h), since this leads to a divergent speed of sound that violates causality. 
One can qualitatively describe a solution as follows; 

a) A solution starts from an initial value h , then, develops towards either a future fixed point 
or a point where F(h) —> ±00. 

b) If it develops towards a fixed point, then, the time taken by a solution to reach this point is 
always infinite according to the argument in section 2 and the solution has no future finite-time 
singularities. 

c) If it develops towards a point where F(h) — > ±00, and the asymptotic behavior is linear or 
slower, then, the time to reach this point is infinite and the solution has no future finite-time 
singularities. 

d) If it develops towards a point where F(h) —> ±00, and the asymptotic behavior is growing 
faster than a linear function but F{h) has the asymptotic form F(h) ~ g/g', where g{h) is any 
function such that g(h) — > 00 as h — > ±00, then, the time to reach this point is infinite and the 
solution has no future finite-time singularities. 

e) If it develops towards a point where F(h) — > ±00, and the asymptotic behavior is growing 
faster than a linear function and F(h) has the asymptotic form F(h) ~ g/g', but g(h) does not 
grow as g(h) — > 00 when h — > ±00, then, the time to reach this point is finite and the solution 
is singular. 

f) In all the above cases the solution is a monotonic function of time as it develops from h 
towards either a fixed point or a point where F(h) — > ±00. 

The same analysis can be followed backward in time but with past fixed points and points 
where F(h) — > ±00. The solution has no finite time singularities if it has no future and no past 
time singularities. 

4 Examples and Application 

In this section we use the qualitative method developed in the previous sections to describe the 
general behavior of viscous fluids in FLRW models and compare it with exact solutions. The 
last exact solution with r = 1/4 describes either as a nonsingular phantom fluid or a unified 
dark fluid. At the end of this section we use the phase space method, causality and stability 
constraints to list possible future scenarios of the universe. 
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4.1 Examples: Nonsingular viscous fluids 

Consider the following equation of state (EoS) 

p(H) = ( 1 -l)p-3r ] (p)H, (18) 

where r)(p) = i] Q p r . The above EoS can describe a fluid with bulk viscosity r)(p) (see e.g., [33]), a 
polytropic fluid [37], or a fluid with adiabatic particle production (see e.g., [36]). Although these 
different interpretations of the above EoS produce the same dynamics, their thermodynamics 
could be different. Several solutions for the above EoS with different values of r, including the 
cases discussed here, are known in the literature, see for example [35]. Here we are going to 
discuss two viscous solutions, the one with r = 1 which is well known in the literature [34] and 
another with r = 1/4, which is less known and express them in terms of Lambert-W function. 
This clearly shows how the density and the scale factor behave as functions of time. 
Now using the above pressure in Eqn. 3 we obtain 

H = -\ (p + p) = ~ (7 H 2 - 3 r Vo H 2r+1 ) (19) 

Taking h = H/H* and r = cH*t, where H* = (3%)^ and c = 37/2 we get 

^ = -h 2 (l-h 2r ~ 1 ), c-^ = h, (20) 
get a ar 

Notice that we have two fixed points, h 1>2 = 0, 1, for this equation. The nature of these fixed 
points depends on the value of r. Solutions with r > 1/2 behave differently compared to 
those with r < 1/2. In the first case F(h) is negative during its interpolation between the 
two fixed points, while the reverse is true for the second case. It is interesting to notice that 
the first case describes a nonsingular universe which has an EoS parameter w(p) > —1 since 
F(h) = — 1/2 (p + p) while the second case describes a nonsingular universe dominated by a 
phantom component with EoS parameter w(p) < —1. It is interesting to notice that fluids with 
r < 1/2 and h > 1 have a similar behavior to that of a generalized Chaplygin-gas cosmology, 
i.e., for large scale factor, it behaves as a cosmological constant and for small scale factor it 
behaves as a fluid with EoS p — (7 — I) p. One can see that by writing Eqn. (4) as 

a ^ = -3 (p(a) + p(a)) = -3 7 p (1 - rf p r " 1/2 ), (21) 

where r]' = \/3 770 solving the above equation for a general viscous fluid with r — 1/2 = — s < 0, 
we get 

p(a) = ( 7 A/r 1/s \l + C ia -M 1,S (22) 
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Figure 3: Fluids with r > 1/2 and r < 1/2 



If p > (V/7) 1 ^ = P*i the integration constant C\ is positive, sice 

Ci = a*' [(p /p*) s - 1] , (23) 
then for small a, i.e., (a /a) 37S [(p /p*) s — 1] >> 1, we get 

p ~ a' 3 '' 

which describe a fluid with an EoS p — (7 — 1) p, and for large a we get 

P ~ P* (25) 

which describes an empty space with a cosmological constant p*. Another interesting feature 
of these models is that they show that a normal matter, i.e., 7 > 0, with bulk viscosity, behaves 
as a phantom matter. Furthermore, in the r < 1/2 case these solutions are nonsingular. To 
show these features, let us start with the pressure 

p = 1P (l-(H*/Hy)-p, (26) 

which leads to an effective EoS parameter 

w eff = 7(1- (H*/HY) - 1, (27) 

since 7 > 0, and H*/H > 1, we always have w e ff < —1, which breaks all energy conditions. 
Since H(t), in this case, interpolates between two fixed points (H* and 0) and the pressure 
p(H) is continuous and differentiable, then according to the argument in section 2, the solution 
is nonsingular and takes an infinite time to reach a fixed point. We are going to see an explicit 
example of this behavior for the r = 1/4 case. 
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(24) 



Fluids with r=l : This case was first discussed in [34] as nonsingular solution for a viscous 
fluid cosmology. In literature this solution usually expressed in terms of time as a function of 
the scale factor or the density. Here we express the energy density and the scale factor as 
functions of time in terms of Lambert W-function. First, let us analyze this case qualitatively 
taking subsection 3.2 in consideration. The asymptotic behavior has the form F(h) ~ h 3 which 
leads to singular solutions unless there are fixed points. F(h) has two fixed points hi — and 
h 2 — 1, as shown in Figure (4). The first is half-stable point and the second is an unstable 
point. These points divide possible solutions into three types; i) a solution where h G (— oo,0], 
if ho < hi, ii) a solution where h G [0, 1], if h 2 > h > hi, and iii) a solution where h G [1, oo), 
if h 2 < h - Notice that, in case b), since p(h) is differentiable then by the argument in section 
2, it takes the solution an infinite time to reach hi starting from some initial value, ho, where 
hi < ho < hi. The same is true if we calculate the time taken by the solution to go from 
hi to ho. Therefore, for hi < h < hi the solution is nonsingular and interpolates smoothly 
between hi, and h 2 . If h < hi or h > h 2 the solution has a finite time singularity since 
the asymptotic behavior of F(h) ~ h 3 is growing faster than a linear behavior. For h > h 2 
the solution describes a universe that starts from a de Sitter space and endes with a Big Rip 
singularity. The above equations can be solved exactly in terms of Lambert W-function, which 



is the solution of the equation W e = x. The hubble parameter and the scale factor in terms 
of the time "t" are given by 



Fixed Points, r=l 




Figure 4: Fixed points for r = 1 case 



H(t) = H*[W(cie^^ HU ) + l}-\ 
a(t) = c 2 [W(cie^ HU )] 2 ^ 



(28) 




(29) 
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a(t) = ao 



W(0) 



(2/3 7 ) 



where (3 is given by 



Mo 



(— -i) 



(30) 



(31) 



It is crucial at this point to know the sign of (3 since it controls the behavior of the W-function. 
The ho = Hq/H* < 1 initial value corresponds to a positive (3, which leads to a smooth 
behavior in all times for the energy density. For the scale factor, the moments at which it 
either diverges or goes to zero are when t = +oo and t = — oo respectively. Therefore, we have 
no finite-time singularities in this case. In early times, this model describes an empty universe 
with a cosmological constant A ~ H* 2 , which evolves to a universe with an EoS p — (7 — 1) p 
in late times. The behavior of H(t) and a(t) as a function of time t, is shown in Figure(5), 
which is clearly monotonic. If the initial value ho = H Q /H* > 1 (/3 < in the case), then we 



ll/thhic /'araiueicr 



1 




o.s\ 




0.6 

HO) 

0.4 




0.2- 






Figure 5: Hubble parameter H(t), in units of H*, and scale factor a(t), in units of a , versus time t in units 
of l/H*, r = 1 and 7= 1. 

have a singularity of type-I, or a Big Rip singularity. In this case the effective EoS parameter 
w e ff = p/p < —1, therefore the fluid is phantom. Notice that, if h < h 2) the r < 1/2 cases 
describe nonsingular phantom solutions. 

One can modify the above model to accommodate the late-time acceleration by adding a 
cosmological constant A. The exact solution of this modified model is not easy to get, but if 
A/if* 2 = A << 1, or the cosmological constant of the early times is larger than that of late 
times, then the fixed points become 



hi ~ VX h 2 ~ 1 — A, 



(32) 



This model describes a universe that interpolates between two de Sitter spaces one with large 
cosmological constant in early times and another with small cosmological constant in late times 
which can model the inflation and late time acceleration periods. 
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Fluid with r=l/4: Here we discuss the solution for the r = 1/4 case, which can be expressed 
in terms of Lambert W- function. As we mentioned above, cases with r < 1/2 are interesting 
since they have a unified dark fluid behavior when h > h 2 . If h < h 2 , the solution describes 
a nonsingular phantom matter with w e ff < —1. To analyze this case qualitatively, let us start 
with the asymptotic behavior of F(h) which has the form F(h) ~ h 2 . It clearly leads to singular 
solutions unless we have fixed points. F(h) has two fixed points hi — and h 2 — 1, as shown 
in Figure (8). The first is an unstable point and the second is a stable point. These points 
divide possible solutions into two types; i) a solution where h G [0,1], if h 2 > h > h±, and 
ii) a solution where h G [l,oo), if h 2 < h . Notice that, in the first case, i), since p(h) is 
differentiable then by the argument in section 2, it takes the solution an infinite time to reach 
hi starting from some initial value, ho, where hi < h < h 2 . The same is true if we calculate 
the time to go from h 2 to ho, therefore, the solution is nonsingular and interpolates smoothly 
between hi, and h 2 . The second case ii) is singular, since the asymptotic behavior F(h) ~ h 2 
is growing faster than a linear behavior. In fact, this solution describes a universe that starts 
from a finite-time singularity in the past (Big Bang type) and evolves to a de Sitter space after 
an infinite time. 



Fixed Points, r=l/4 




Figure 6: Fixed points for r = 1/4 case 



Eqn.(19) for r = 1/4 case has an exact solution which has the following form 



H(t) 



H* [W(ci e 



(-3 7 /4) H*t 



) + l] 



-2 



a(t) 




(33) 



Using initial conditions, the integration constants are 



4 





(34) 
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where 



(35) 



Notice that f3' is positive for /i < 1 and negative for /i > 1- F° r initial value h < 1 one expects 
a nonsingular solution which interpolates monotonically between Minkowski space and de Sitter 
space. The Hubble parameter H(t) and the scale factor a(t), as functions of time "t" , are shown 
in Figure(7). As we mentioned in the beginning of the section, this class of solutions shows how 




Figure 7: Hubble parameter H(t), in units of H*, and scale factor a(t), in units of a , versus time t in units 
of l/H*, r = 1/4 and 7 = 1. 

a viscous fluid with a usual EoS (i.e., 7 > 0) behaves as phantom component. Similar to the 
r = 1, one can consider adding a small cosmological constant. A small cosmological constant 
A/if* 2 = A << 1 leads to the following fixed points 



/ii~A 2/3 , /i 2 ~l-2A. 



(36) 



The solution with initial value hi < h < /i 2 described a universe filled with a phantom matter 
interpolating between a de Sitter space with a small cosmological constant at early times and 
another with a large cosmological constant in late times. 



4.2 Fate of the Universe 

Here we list all possible future scenarios of the universe as a single fluid in FLRW cosmology 
without assuming any fixed points but imposing the causality and stability constraints. It is 
known that, if p > —l/3p, destiny of the universe is tied to geometry (see for example [38]) and 
the value of k is important to predict the fate of the universe. On the other hand, iip < — 1/3 p 
destiny is not tied to geometry but controlled by the behavior of the energy density p. This can 
be shown if we consider the known mechanical model for a, by rewriting Friendmann equation 
as 



d 2 



a 2 p(a)-k = -2V eff { 



(37) 
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Taking the EoS p = w p leads to 

p(a) = Ca- 3(1+w) =>2V eff (a) = -C /3 a~ (1+3w) + k (38) 

For large a, if w > —1/3, k controls the existence of vanishing velocities, or turning points, 
but for w > —1/3, the potential gets a small contribution from k compared to that coming 
from a 2 p. This breaks the connection between the geometry and destiny. In fact, this simple 
argument is also suggesting that our universe will keep on expanding because of the domination 
of the dark energy component. We will see next that this is not generally the case. 
Here we are going to use the above constraints to list possible scenarios for the future of the 
universe. Let us model our universe using a general single barotropic fluid, which is a reasonable 
assumption since dark energy is dominating in late times. 

H = 1/2(3 H 2 +p(H)). (39) 

First, it is easy to show that we are in a region in the phase space (i.e., H — H space ) where 
H < 6 . It is known that the deceleration parameter has changed sign from positive to negative 
as the universe evolved from a matter dominating era to a dark energy dominating era. To see 
how this crossing happened consider first the zero acceleration curve which is given by 

- = H + H 2 = H = -H 2 , (40) 

a 

then, as an example of dark energy, consider an accelerating universe with a cosmological 
constant and matter 

H = F(H) = -l/2(3H 2 -p A ). (41) 

By plotting both functions — H 2 and —1/2 (3H 2 — p A ) in Figure(8) one can see how this crossing 
happened, i.e., going from a < to a > 0. The point at which the acceleration vanishes can be 
used as a reference point to draw the constraints on possible evolutions of the universe. Using 
the causality and stability constraints, dp/ dp < 1 and dp/ dp > 0, we get > dp(H)/dH <6H 
which in turn leads to —3H > dF(H)/dH < —6H. Integrating this inequality leads to 
—3/2H 2 + C\ > F(H) < —3H 2 + C2, where the integration constants C\ and C2 are fixed by 
the initial values of H and H at the reference point. Notice that, C\ and C 2 have to be positive 
numbers, otherwise crossing the zero-acceleration curve will not occur. The last inequality 
constrains the future behavior of F(H) to lie between these two parabolas, as a result, F(H) 
must meet the if-axes in a future time and ends as a de Sitter universe after infinite time. 
In addition, this evolution starting from the point of zero-acceleration till the end of time is 
nonsingular since it ends with de Sitter universes. 

6 One can directly calculate H using the best measured values for w Dl TO and Q, D . 
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Figure 8: The —H 2 curve is the zero-acceleration curve. The black region between the two curves, —3/2H 2 
and —3/2H 2 satisfies the two constraints on barotropic fluids and any curve in this region must end with a 
fixed point on positive iJ-axes 

As one might notice, although, the causality constraint is essential for any physical model. 
It is not clear if we should insist on having the stability constraint, since we do not know the 
physics of the dark energy component. Now, if we relax the stability constraint, dp/ dp > 
and allow the pressure derivative to go negative. From the above phase space diagram, it is 
clear that the universe either ends as an empty universe or hits the negative if-axes. The 
last possibility is interesting since according to the discussion in subsection 3.1 it describes a 
turnaround behavior, therefore, the universe in a future finite time reaches a maximum size, 
then, recollapses, since the hubble parameter H changes sign. 

5 Conclusion 

In this work, we used a phase space method to study possible consequences of having fixed 
points in a single fluid flat FLRW models. Some of these are; (i) if we describe our universe as 
a single component fluid with a future fixed point, then the resulting cosmology does not have 
future-time singularities of types I, II and III in [25], (ii) cosmologies with a future and a past 
fixed points are free of of types I, II and III singularities, (in) one can use a simple argument 
to show the phantom divide [27, 28], or in a single fluid FLRW models it is impossible for a 
physical solution to cross the phantom divide line in a finite time, and (iv) in these models, the 
only way to get bounce solutions is to have a nonvanishing pressure as p — > 0. This method can 
be used to construct nonsingular late-time models, in particular, unified dark fluid and dark 
energy models. We use this method to qualitatively describe any flat FLRW model with fixed 
points. We discussed FLRW cosmology with bulk viscosity rj ~ p r , and presented two exact 
solutions with r = 1 and r = 1/4, which are expressed in terms of Lambert-W function. The 
last solution describes either a nonsingular phantom dark energy or a unified dark fluid model. 



17 



The phantom solution is interesting since it shows how a viscous normal fluid behaves very 
similar to a phantom matter without Big Rip singularities. In addition, it interpolates between 
two de Sitter spaces with small and large cosmological constants. Possible future scenarios of 
our universe include; a de Sitter space, an empty universe with vanishing cosmological constant, 
or a turn a round solution that reaches a maximum size, then collapses. 
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